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Abstract 

At short distances, energy eigenfunctions of chaotic systems have spatial 
correlations that are well described by assuming a microcanonical density 
in phase space for the corresponding Wigner function. However, this 
is not correct on large scales. The correct correlation function is in 
turn needed to get the correct formula for the root-mean-square value 
of the off-diagonal matrix elements of simple observables, and for the 
fluctuations in the diagonal elements. 
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My focus in this talk will be on the statistical properties of the quantum energy eigen- 
functions in classically chaotic systems. Twenty years ago, Berry conjectured that these 
eigenfunctions could be treated as gaussian random variables, and that the spatial corre- 
lations would be those following from taking the expected value of the Wigner density for 
each eigenfunction to be microcanonical. More specifically, the probability that the actual 
eigenfunction V'o(q) for energy E = Ea (the corresponding energy eigenvalue) is between 
ip{q) and ip{q) + dtplq) for all coordinate points q is given by 



P^iplE) oc exp 
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. (1) 



If the system is time-reversal invariant, the eigenfunctions are real, (3=1, and the measure 
is "DV = Hqdip^q); if it is not, the eigenfunctions are complex, /3 = 2, and the measure is 
Vip = Y\^dlle'ijj{q) dlmiplq). In either case, the kernel K[q^,c['\E) is the functional inverse 
of the two-point correlation function 

C(q',q|E) = (^(q')^*(q)), (2) 

where the angle brackets denote averaging over P{iIj\E). Taking the Wigner density to be 
microcanonical results in 

^'^^ = / (2^ e^^-^^'-^^^'SiE - Hwip, q)) , (3) 

where q = |(q+q') is the midpoint, H\y{p, q) is the classical hamiltonian (more specifically, 
it is the Weyl symbol of the hamiltonian operator), and p{E) is the semiclassical density of 
states, 

P{E) = Jj^^(E-Hw{p,ci)). (4) 

Berry's conjecture for C(q', ql-E) can be tested numerically in two-dimensional billiard 
systems; in this case, H = p^/2m, and we have C{q_',q\E) = Jo{k\(l' — q|), where E = 
fi^k'^ /2m and Jo{x) is an ordinary Bessel function. An "experimental" correlation function 
Cexp(s, D\E) can be defined in terms of a particular energy eigenfunction ^/'^(q) via 

Ce.p(s, D\E)^^ I d\ rM + - ¥) , (5) 

where E = E^, As is the area of the billiard, and is the area of a domain D over 
which the central point q is averaged. Early computations of Cexp(s, -D|-E) found that 
it typically resembled its predicted value Jo(fc|s|), but with large fluctuations. Later it 
was realized that these fluctuations were in fact entirely consistent with the expected 
gaussian fluctuations ^ of V'o(ci)- Recent detailed computations leave little doubt 
that (at asymptotically high energies) the eigenfunctions of chaotic billiards are very well 
described by Berry's conjecture. 

Berry's conjecture also follows from other approaches to quantum chaos. The random- 
matrix analogy P,^] suggests that the energy eigenstates have a Porter-Thomas distribution 
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of amplitudes in a suitably chosen basis. For billiards, the natural choice is the momentum 
basis, and from this starting point one can derive eqs. (|l]) and (^) . If one adds a white- 
noise, spatially uncorrelated random potential to the system, and takes the limit where the 
mean-free path is much larger than the system size, then eq. (|l|) can be derived via the 
supersymmetric sigma-model technique [^1|-|16|. 



However, we cannot expect eq. (y) to be valid when the separation between q and q' 
becomes large compared to distance scales which are important classically. This is because 
the hamiltonian evaluated at just the midpoint q does not contain enough information about 
the classical landscape between q and q'. To find the correct formula for C(q',q|i?) when 
|q' — q| is large |T^, we first consider the theory with an added random potential. In this 
case an expression for C{q^,q\E) can be derived explicitly [|T2| : 



C(q', q|^) = [^(q, c^lEY - G(q', q|E)] . (6) 

Here G(q',q|-E') is the energy Green's function, 

G(q,q|i5)-L^_^^+,0+' 
p{E) is the density of states, 

p{E) = ^Jd^Q [^(q, q|i?)* - G(q, q|E)] , (8) 

and the bar stands for averaging over the random potential. We now make the assumption 
that, for a chaotic system, eq. (|^) continues to hold if we first take the limit of small h, and 
then disregard the random potential. 

There are two different semiclassical approximations which can be used to compute the 



Green's function in the small-^ limit [p!8| , |l9| . The first applies when the distance between 
q and q' is small, in the sense that the classical path of least action with energy E which 
connects q to q' is well approximated by a linear function of time (in any coordinate system). 
This path then dominates the Green's function, and one finds 

G(q q\E)= f ip{q'-q)A I ('9^ 

which immediately yields eq. @ for C(q', ql-E"). If, however, the classical path of least action 
is not well approximated by a linear function of time, and the value of this action is much 
larger than h, then we have instead 

11^' - M2Jn)u-m S W^e'^'"-'"' . (10) 

Here the sum is over all classical paths connecting q to q' with energy E, action 

Sp = r p-dq, (11) 



q 
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focal point number z/p, and fluctuation determinant 

(12) 



Dp = det 



d'^Sp d^Sp 
\ d(i2dE 



Eqs. d^) and (|T0|) give the correct formula for C(q', c\\E) when |q' — q| is large JIT . 

This formula is most useful under circumstances where the sum over paths is dominated 
by a single path of least action. An interesting example is a two-dimensional billiard in a 
perpendicular magnetic field B. In the billiard interior the hamiltonian is H = (p— eA)^/2m, 
and we will work in the gauge in which the vector potential is A = |B x q. Assuming 
that it is not blocked, the classical path of least action is a circular arc with length i, 
related to the separation L = |q' — q| and classical cyclotron radius R = {2mEY^'^ /\eB\ via 
£ = 2i?sin~^(L/2i?). The action for this path can be divided into a geometric part and a 
gauge-dependent part, Sp = Sgeom + •S'gaugo- The geometric part is 

'S'gcom ^ 

where hk = {2mEf/'^, and A=\m- ^R^sm{i/R) is the area enclosed by the circular arc 
and the straight line connecting q to q'. For our gauge choice, the gauge- dependent part is 

-^gauge = |eB-(q X q') . (14) 

In any gauge, the determinant \Dp\ is given by 

Keeping only the contribution of this one path, and recalling that p[E) = mAB/27ch'^ for a 
two-dimensional billiard with area As, we find from eqs. and (|1^) that 

We need kL ^ 1 for this formula to hold. Also, to avoid an integrable region of phase space, 
a circle of radius R must not be able to fit inside the billiard; this implies L < 2R. 

We can compare eq. ([T6|) with the result of using Berry's formula, eq. (^; with our gauge 



choice we find 

C(q',q|E) = A^'exp{iSg^^gjn)MkL) . (17) 
For large kL, we can use the asymptotic form of the Bessel function to get 

C(q',q|E) = ^B'exp(z5„//^)^^|^y^ • (18) 
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We see that Berry's formula misses the corrections due to the finite cyclotron radius R which 
are present in eq. ([16|) . 

If we make a gauge transformation 

A(q) ^ A(q) + V<l'(q) , (19) 

where $(q) is any smooth function, then 

5'gaugc ^ 5'gaugc + e$(q') - e<l>(q) , (20) 

and so eq. ([I6|) implies 

C(q', q|E) ^ e+*^[*('i')-*('i)l/^C(q', q|E) . (21) 
That this is correct can be seen by recalling that a wave function ip{q) transforms as 

^(q) ^ e+'^^^'i^/Vlq) (22) 



under (|T9|), and that C(q',q|£'a) is the expected value of '?/'Q(q')V'a(q)- the other hand, 
eq. (^ implies instead that 

C(q', q|E) ^ e+''^^'-^y^'fin)/nc(^^^ q|^) ^ (23) 

which again illustrates the fact that Berry's formula is valid only when |q' — q| is sufficiently 
small. 

I now turn to another issue, the statistical properties of transition matrix elements Aa/B = 
{a\A\i3) in chaotic systems. Here the operator A is a smooth, ^-independent function of the 
coordinates and momenta, and |a) and are different energy eigenstates. We will be 
interested in the average value of |^o/3p when the energies Ea and are each varied over 
a range which is small classically, but encompasses many quantum energy levels. For later 
convenience let us define 

E = ^{Eo, + Ep) and huj = Ep-E^. (24) 

We will do our computations in the limit % ^ {] with E and uj held fixed. The mean level 
spacing near energy E is l/p{E) ~ , so our condition on the ranges of Ea and Ef^ is easy 
to fulfill. 

One approach jS^-^ to this problem is to relate the average value to the 

classical time-correlation function of the observable A via 



Aap\^) = — r dte''^' f d^^EAw{Pt,(lt)Aw{p,(l) ■ (25) 



Here th = 27ihp{E) is the Heisenberg time, (pt, q^) is the point in phase space which is 
reached classically at time t when starting from (p, q) at time zero, Aw[p,q) is the Weyl 
symbol of the operator A, and dfiE is the Liouville measure on the surface in phase space 
with energy E, given by 
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Note that the notation, though useful later, is somewhat misleading: djj^Eisa purely classical, 
^-independent object. We will not discuss the derivation pO|-^ of eq. (|^) here; a brief 



and non-rigorous review is given elsewhere ^3 

Another approach p|,^|7| to computing (lA^^p) is to write in terms of the energy 

eigenfunctions V'alQ) ^ind ipf3{ci), and then average |Aa/3p over the eigenfunction probability 
distribution of eq. (|lD. This calculation is simplified if A is a function of q only (rather than 
both q and p), and so we specialize to this case. We get 



(27) 



It is not at all clear that eq. (^) gives the same result for (IAq,^! ) as eq. (pSD, an issue which 
was raised (in the context of the diagonal matrix elements) by Austin and Wilkinson p3 . 
If Berry's formula is used for C(q',q|-E), then eqs. (^) and (pTf ) do not necessarily agree, 
as can be seen by working out some simple examples. In this case, it is eq. ( [27| ) which is 
wrong. The reason is that q and q' are independently integrated in eq. (pT]) , and so |q' — q| 
is generically large. Therefore we should use eqs. @ and (|10D for C(q',q|i?), rather than 
eq. (H). If we do, then eqs. (p5|) and (|27| ) give the same result for (|Aq,^P). 

To see this still requires some work. A detailed exposition is given elsewhere [23], and 
here we will just highlight the key elements. One is the diagonal approximation after 
substituting eqs. @ and ([T0|) into eq. (|27D , the double sum over paths is collapsed to a 
single sum, since the off-diagonal terms will have rapidly oscillating phases. In different but 
related contexts p6| , p^ , this requires restricting the single sum to paths whose elapsed times 
are less than the Heisenberg time, and we will assume the same is true here. The action 
difference in the single sum is 



huJTr, 



E. 



dE 



+ ... 



E=E 



(28) 



where Xp is the elapsed time along the path. We assume that i/p^ = i/p^, 
topological quantity which in general will not change when the energy of the path is varied 
slightly. Finally, we note that the determinant Dp can be written as 



= det 



9p 


dr 








dr 


dE 


dE 



(29) 



Here p = ~dSp/dq is the momentum at the beginning of the path, and r = Tp is the elapsed 



time along the path. Eq. (^) shows us that \Dp\ can be thought of []TB[ as a jacobian for a 
change of variables from the final position q' and total energy E to the initial momentum p 
and elapsed time r. Eq. (^) already has an integral over q', and to get one over E we insert 
1 = / dE 6{E — Hw{p, q))- Now we can make the change of integration variables suggested 
by eq. (|29|) , and we find 



Aap\') 



1 



nhp{Ey 



dr 



d^p d^q 



^ 5iE - Hwip, q)) cosiiUT)Aw{q!)Aw{ci) . (30) 



paths 
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The sum is now over all paths which begin at (p, q) and have elapsed time r. However, 
there is only one such path, and so the sum over paths may be dropped. Also, q' is the 
position at time r, and it is now more properly denoted q,-. Using eq. (pB]), the fact that 
time-translation invariance implies that / dfi^ A.w{(It)A.w{^) is an even function of r (even 
if the system is not time-reversal invariant), and th = 2nhp{E), we see that eq. ( pOD can be 
rewritten as 



1 Z'+'^H 



(31) 



which is equivalent to eq. (|25|). 

Another quantity of interest is the size of the fluctuations in the diagonal matrix elements 
Aaa- If we first shift A (if necessary) so that {Aaa) = I dfiE^ AvK(q) = 0, then the object 
we wish to evaluate is This has been done previously by making use of 

the trace formula |P^pU| , pU[ and properties of periodic orbits. Here we will instead compute 
(|^Q«P) by averaging over the probability distribution for energy eigenf unctions P,p4|,|7|,p3 



In the case of a system which is not invariant under time reversal, the energy eigenfunctions 



are generically complex, and the relevant formula is [15 



(^l*7A2^3>4) = (^3>4) + (^3*^2) , (32) 

where i/ji = il^a{(li)- If the system is invariant under time reversal, the energy eigenfunctions 



are real, and we have instead |T5| 

(^1^2^3^4) = (^"3^4) + (V'lV'4)(V'2^3) + (V'lV'3)(^2V'4) 

Combined with the previous results for (IAq/jP), we find that 



2/P 



dr j dpE^Aw{(lT)Aw{(l) ■ 



(33) 



(34) 



Here /3 = 1 for a system which is invariant under time reversal, and {3 = 2 for a system 



which is not. Eq. ([3^ ) is in agreement with the earlier results p7|j28|j24 . 
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